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Abstract. We analyse theoretically for the first time to our knowledge the perfect 
phase matching of guided TE and TM modes with a multilayer waveguide composed 
of linear isotropic dielectric materials. Alongside strict investigation into dispersion 
relations for multilayer systems, we give an explicit qualitative explanation for 
the phenomenon of mode matching on the basis of the standard one-dimensional 
homogenization technique, and discuss the minimum number of layers and the 
refractive index profile for the proposed device scheme. Direct applications of the 
scheme include polarization-insensitive, intermodal dispersion-free planar propagation, 
efficient fibre-to-planar waveguide coupling and, potentially, mode filtering. 

As a self-sufficient result, we present compact analytical expressions for the mode 
dispersion in a finite, iV-period, three-layer dielectric superlattice. 
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1. Introduction 

Since the early experiments on mode conversion in film-waveguide magneto-optical 
systems [H |2] and the subsequent theoretical works [3 IH El it has become clear that an 
equality between propagation constants of the two guided modes — otherwise called mode 
matching or perfect phase matching — can be realized exclusively due to the tensorial 
character of dielectric permittivity of the waveguiding slab, i.e. due to its either natural 
or induced anisotropy. An obvious way to gain anisotropy is to produce a waveguide 
layer of mono crystalline material possessing intrinsic birefringence; another possibility 
is provided by use of gyrotropic waveguides (for example, YIG on GGG) subject to an 
external DC magnetic field. Historically, it was the latter way that experimentalists 
followed, probably because of its higher flexibility and technical reliability. A number 
of device schemes implementing that principle were rapidly proposed, from an optical 
switch to a sensitive experimental method of determining the permittivity tensor 
elements The two unpleasant drawbacks were then the complexity of the schemes 
and the comparatively high absorption in gyrotropic films that hindered their utilization, 
e.g. for intermodal dispersion-free planar propagation. 

On the other hand, in recent years a certain interest in traditional multilayer 
dielectric mirrors from the viewpoint of photonic bandgap localization has been 
observed, stimulated probably by the unexpected discovery of total omnidirectional 
reflection from an appropriately designed one-dimensional dielectric lattice, i.e. an 
actual photonic band gap exhibited in case of such a simple geometry [HI El Ej ■ What 
looks now curious but is important for the purpose of the present treatment is that one 
can easily see the points of perfect phase matching of TE and TM modes on the standard 
band structure diagrams of dielectric multilayers [lOj ch 4], yet unfortunately little or 
no attention has been paid in the literature to those mode matching phenomena — in 
spite of so many promising applications including polarization-insensitive, intermodal 
dispersion-free planar waveguiding, efficient fibre-to-planar waveguide coupling and, 
potentially, mode filtering, drastic improvement of mode conversion efficiency and 
diffraction efficiency in magneto-optical Bragg cells and others. To the best of our 
knowledge, the underlying idea that a multilayer waveguide can be regarded as a 
homogeneous, effectively anisotropic slab — and as if it were the case of a naturally 
birefringent or gyrotropic, magnetized waveguide, the perfect phase matching of TE 
and TM modes must occur — has been formulated just recently [TTj . 

The aim of the current paper is to highlight and theoreticalUy investigate the perfect 
phase matching of guided TE and TM modes with a multilayer device composed of 
linear, isotropic, non-magnetic and non-absorbing materials. The paper is organized 
as follows. ISection 2l is devoted mainly to qualitative models and considerations: we 
begin with a transparent explanation for the mode matching phenomenon on the basis 
of conventional one-dimensional homogenization method and then discuss the general 
eigenproblem for TE and TM modes in order to derive some basic requirements for the 
scheme, such as the minimum number of layers and the proper refractive index profile. 
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Strict analytical expressions for the modal dispersion structure of periodic dielectric 
stacks with two-layer and three-layer unit cells follow those qualitative speculations and 
are given in lsection 31 Finally, we analyse in lsection 41 several examples of perfect phase 
matching with simple planar structures composed of silicon and silicon-nitride layers, 
fully compatible with modern integrated optics techniques. 



2. What is mode matching by a multilayer dielectric waveguide? 

A straightforward way to make the idea of perfect phase matching transparent is to 
consider a stack of isotropic dielectric layers as an effectively homogeneous, uniaxially 
birefringent medium, by means of the standard one-dimensional homogenization 
technique. This approach allows to expand the argumentation for the mode matching 
with anisotropic waveguides to the case of dielectric multilayers, with an important 
consequence that the mode matching is in principle possible with those structures too. 
To find the minimum number of layers required for the matching and to optimize the 
refractive index profile, we furthermore analyse the waveguiding conditions for TE and 
TM modes in multilayer structures, tracking the fruitful analogy with confined one- 
dimensional electron states problem quite profoundly investigated in the framework of 
classical quantum mechanics and solid-state physics. 



2.1. Effective medium treatment 

In a single-layer isotropic waveguide, the two given modes — TE (s-polarized) and TM 
(p-polarized) — are commonly characterized at a fixed frequency by definitely different 
propagation angles 6s,p between the wavevector and the plane of the guide, and 
consequently, by different (non-equal) propagation constants I3s,p via 

I3s,p = n cos 9 s^p, (1) 

where n is the refractive index of the waveguiding medium and the subscripts s and 
p differentiate between TE and TM modes. Anisotropy of the waveguide accounts 
evidently for the following modification of equation (1): 

I3s,p = ns,p{0) cos 9s,p, (2) 

where the indices 'seen' by TE and TM modes, and a situation like that in 

figure 1| becomes possible, when the projections of the wave vectors of different modes 



onto the plane of guidance coincide, giving rise to the phenomenon of mode matching. 

Consider now a two-component periodic multilayer system consisting of layers of 
dielectric permittivity ei = n\ and thickness ti separated by layers with permittivity 
€2 = n\ and thickness ^2- If a plane wave is incident onto this array with its electric vector 
E polarized parallel to the layers (TE-wave) and the thicknesses ti and t2 are small 
compared to the wavelength, then one may introduce the effective dielectric permittivity 
of the array as [121 

ti ei + ^2 ^2 /o\ 



Dielectric multilayer waveguides for TE and TM mode matching 



4 




and if the magnetic vector H is parallel to the plates (TM-wave), then the effective 
permittivity becomes ^21 

= ^AlMim. (4) 

ti Ci + ^2 ^2 

This treatment allows one to state that two isonormal, s- and p-polarized waves 
'see' in general different refractive indices of effective waveguiding medium, hence if to 
start from the isotropic (ni = ^2) case and then change one of the indices, the families of 
TE and TM dispersion curves will undergo non-equal shifts on the dispersion structure 
diagram and thus intersect at a number of points, displaying the modes degeneration — 
the perfect phase matching. One should remember by the way that it is of course not 
complete, spatial degeneration of modes, but mere propagation constant coincidence. 

The question that naturally arises here is: How much layers are enough for the 
mode matching to occur? In other words: When does the reasoning based on the form 
birefringence approach lose its validity? To answer these questions and to elucidate 
the physics of mode matching in the case of a finite number of layers, one may find 
it convenient to employ an explicit mathematical analogy between confined photon 
and electron states, the latter being quite thoroughly investigated in classical quantum 
mechanics and solid-state physics — a comprehensible ID Kronig-Penney model should 
certainly be mentioned here as an example. 



2.2. Dispersion equation for TE and TM modes as an eigenproblem: Qualitative 
analysis 

As is commonly known and excellently described by Joannopoulos et al fOl ch 3], 
symmetry considerations allow us to separate the modes of planar dielectric waveguides 
into two classes: TE and TM polarized. In a coordinate system in which the x-axis is 
normal to the bimedium interfaces and the light propagates along the 2;-direction, we 
can write for TE and TM polarizations, respectively, 

E{r) = exp{i-(3z)(f)^{x)y (5) 
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and 

H{r) = exp(i5^/?z)0^(x)y, (6) 

where ui is the angular frequency, c is the vacuum speed of hght, (3 is what we 
call the propagation constant, and (j)^ (x) are the unknown field distributions. 

Substituting one of these fields into corresponding Maxwell wave equation, in the case of 
permittivity e(r) = e{x) being a step-wise function equal to n'j for the jth homogeneous 
layer we obtain the following eigenproblem: 

0. (7) 

This is the dispersion equation for a multilayer waveguide; the difference between TE 
and TM polarizations lies exclusively in different boundary conditions at the interfaces 
of the layers. Thus, at the interface x = xq we have: 

Conditions for the derivatives originate from continuity of 2;-components of magnetic 
and electric field vectors via Maxwell equations for V x £^ and V x H respectively. 

Note that the scalar Helmholtz wave equation (7) with boundary conditions for 
the TE mode (8) is isomorphous to the familiar stationary one-dimensional Schrodinger 
equation 

^^ + ^(^-f/(x)Mx) = (10) 

for a particle of mass m in potential energy profile U{x), £ being the energy levels 
corresponding to eigenfunctions 'ipei^x). 

The difference in boundary conditions for the derivatives of (j)^{x) and (j)^{x) 
functions leads to strict alternation of TE and TM dispersion curves of one- and two- 
layer waveguides. In other words, at a fixed frequency u the propagation constants 
of confined TE modes are always separated by those of TM ones. This 'rule' can be 
illustrated by |figure~2 for the lowest modes in a single- layer waveguide: the ratio X/2d {X 
is the wavelength within the guide) for the (f)^{x) function of the TEq mode is obviously 
more than that for the TMq mode because of the steeper slopes of the 0^(x) function 
near the interfaces due to the boundary conditions, but however high the refractive index 
of the guide and consequently the slopes of (j)^{x) function in vicinity of the interfaces 
may be, the ratio \/2d for the TMq mode will never drop below unity, while for the 
next TE mode — TEi — it is definitely less than 1, actually lying between 1/2 and 1. 

Similar argumentation still holds for the modal dispersion structure of a two-layer 
stack, no matter of the refractive index contrast of the layers, and, in general, for A^-layer 
stack with refractive indices of the layers regularly increasing or decreasing throughout 
the array {rii < n2 < ■ ■ ■ < oi rii > n2 > ■ ■ ■ > n^). 
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Figure 2. Illustration to the rule of alternation of TE and TM modes in a single-layer 
waveguide. 

Surprisingly, a three-layer waveguide of step-wise refractive index profile develops 
considerably different from the above case, provided the second layer is of lower refractive 
index then its surroundings: n2 < rii^^. Also assume for the sake of simplicity rii = 
and, as before, symmetric sandwiching: no = 1 for both ambient media. As it will now 
be argued, the mode matching can be effectively implemented with such a system; to 
understand its origin, it might be a helpful trick to recall the changes in the system of 
eigenfunctions ipsi^) of electron in one-dimensional three-step potential gap, each 
step being of width d/?>: 

{Uq ioT \x\ > d/2 (outside the gap) 
U2 for |x| < d/6 (in the middle region) (11) 
otherwise (in the lateral regions) 

with increasing U2 from U2 = t/1,3 to U2 = Uq that corresponds to decreasing 77,2 from 
^2 = ^^1,3 (isotropic single-layer waveguide) to unity (two waveguide channels separated 
by air region), the frequency uj being fixed. 

While the difference between and U2 is negligible, the energy levels 

correspondent to ipslx) eigenfunctions compose a relatively regular discrete energy 
spectrum, and so do the propagation constants of TE modes at any fixed frequency. 
The TM mode spectrum is a bit shifted with respect to the TE one due to the difference 
in boundary conditions for the modes of different types, and as just emphasized, the 
propagation constants of TE and TM modes are strictly alternating. 

By contrast, when U2 is large enough, the regions of low potential become weakly 
coupled, the entire system possessing the twofold (generally A^-fold, where is the 
number of regions of low potential — or high refractive index) splitting of mixed odd 
and even energy levels corresponding to nearly degenerate TE modal groups of A^ 
modes each. TM modes also form similar A^-fold groups which are, however, shifted 
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Figure 3. (a) Schematic diagram of the three-layer dielectric waveguide on silica 
substrate and (b) its refractive index profile. 



with respect to TE groups. It is clear now that for certain n2 lying between those 
extremes, or equivalently, for certain frequency uj in case of 77-2 being somewhat less 
than ni^3 and fixed, the matching of TE and TM modes should be expected, with the 
only exception for TEq mode in agreement with the inevitably negative character of 
form birefringence We can also anticipate herefore that the minimum sufficient 

number of layers for the perfect phase matching is three, provided the central layer is of 
lower refractive index than the surrounding ones — see [figure 3| for an example of such 
a structure. More generally, the refractive index profile of the waveguide should have 
more than one maximum for the effective mode matching. 

Consider now a kind of three-layer system 'reciprocal' to the above — with a middle 
layer being of higher index than the lateral ones. For this third possible arrangement of 
the layers, the mode matching is not forbidden, yet there is seemingly no solid ground to 
expect the matching of the lowest modes now since such a system appears to stand close 
to conventional slab waveguide — the latter with no possibility for the mode matching 
at all. More thorough analysis in the manner above for a single-layer waveguide shows 
that coincidence of propagation constants should be expected for the modes of the same 
order ^ but as before with the exception of TEq and TMq modes. It is demonstrated in 
Isection 4l that the mode matching is indeed observed in this case only in the multimode 
regime, thus actually losing its practical value. 



3. Dispersion analysis for multiperiod two- and three-layer waveguides 

The general form (7)-(9) of the dispersion equations for guided modes is quite good 
for qualitative analysis, but inconvenient for numerical investigations into the TE and 
TM mode spectra of multilayer stacks. In the current section we render the known 
dispersion relations for multiperiod bilayer waveguides and present their generalization 
for three-layer multiperiod systems. Thus, following Yeh Jl], in case of an A^-period 
stack with a bilayer unit cell one has for TE modes: 

[f J_ + J_) cc - f-^ + ss 
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rjj = ^rij — (3'^, C,sub,Cov = y P'^ ~ ^sub,Cov ^"^^ shorthand notation such as 
SC stands for the corresponding harmonic function product, in the case above, 
sin(Hini) cos(fct2n2), with the free-space wave number k = u/c. 
Likewise for TM modes IT 
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(15) 



The above equations look rather awkward, but they can be easily handled and 
promptly solved numerically in any mathematical computing package. 

Three-layer semi-infinite dielectric superlattices have been studied earlier, but with 
the focus on the existence and the dispersion of bulk and surface polaritons In the 
case of a finite, iV-period stack with a trilayer unit cell one can however quite trivially 
derive the following dispersion relation for the guided TE modes (see the appendix! for 
the details): 
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(16) 
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and for TM modes: 
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(19) 

,%ei '7ie3y 2 \ 77362 ?72e3 / 2 
Of course, certain symmetries could obviously be found in these equations, but still 
the latter seem to be too cumbersome for any qualitative speculations. Instead, one can 
easily do now plenty of computational work related to physically realizable multilayers, 
and in the following section we will do that with waveguides made of such popular 
optical materials as silicon and its compounds. 



4. Simulated results for multilayers of silicon and silicon nitride 

To demonstrate the performance of dielectric multilayers for the mode matching and 
to corroborate theoretical considerations presented in Isection A we have plotted in 
[figure 4| the modal dispersion curves of (a) two-layer Si/Si3N4 and (b) three-layer 
Si/Si3N4/Si waveguides, as well as of (c,d) their two-period counterparts. Silicon and its 
nitride have been chosen deliberately since these materials are found to be suitable for 
manifold passive optical components like branching waveguides, couples, interconnects, 
interferometers, power splitters, and filters fOJ El- We employed here the mode 
dispersion relations for bilayer and three-layer waveguides and their generalization for 
multiperiod systems, equations (12)-(19). The thicknesses of the layers were assumed to 
be equal; the refractive indices are 3.5 (Si), 2.0 (Si3N4) and 1.5 (Si02). All the numerical 
calculations were performed with Mathematica 4.1 |18j . 
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Figure 4. TE ( ) and TM ( ) families of tlie modal dispersion curves for (a) 

bilayer Si/Si3N4, (b) three-layer Si/Si3N4/Si, (c) two-period four-layer (Si/Si3N4)^ and 
(d) five-layer (Si/Si3N4/Si)^ waveguides, all on silica substrate. 



Figure 4(a) is of no particular interest — we see the system of TE and TM dispersion 
curves that is quite typical for slab dielectric waveguides, with just slight perturbations 
of the high-order modal curves near (3 ~ 2 — that is close to the refractive index of 
the middle layer. In contrast, figure 4(b) reveals the predicted intersection of modal 
dispersion curves for the three-layer system; for the first phase-matching point where 
the TMq mode matches TEi at NT = 0.25A (here T is the total thickness of one period, 
so that NT gives the thickness of the whole structure), we can estimate the thicknesses 
of the layers to be ti = ^2 = ^3 = 130 nm for operation near 1.55 /im, i.e. in the 
spectral range for telecommunication frequencies. In order to improve the operation 
bandwidth, the slopes of dispersion curves at the given point of intersection must be 
as close as possible, and for that purpose the refractive index distribution should be 
further optimized. 

The comparison between figures 4(b), (c) and (d) confirms that it is the number 
of layers of high refractive index, and not periods — which sometimes are more 'mental' 
constructs than 'physical' objects — that determines the number of modes in the TE and 
TM bands formed. Thus, for both Si/Si3N4/Si and (Si/Si3N4)^ structures constituted of 
two silicon layers — two 'potential gaps' for photons — the twofold grouping of modes at 
high j3 values is clearly observed. Finally, [figure 5| testifies to the statement of Isection^ 
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Figure 5. Modal dispersion curves for three-layer Si3N4/Si/Si3N4 waveguide on silica 



substrate — the 'reverse' structure of that in figure 4(b) 



that for the effective mode matching, the refractive index of the medium layer in three- 
layer waveguide should be low enough; otherwise, the matching is not forbidden, but 
occurs only for high-order modes and remains thus impracticable. 



5. Conclusion 

In summary, the perfect phase matching of TE and TM guided modes with a multilayer 
waveguide composed of linear isotropic dielectric materials has been anticipated and 
simulated for the first time. Dielectric superlattices with unit cells of two types have been 
considered: two-layer and tree-layer; for finite, A^-period, three-component structures, 
the compact analytical expressions for the mode dispersion have been derived, which 
can be regarded self-sufficient result. 

Numerical simulations performed for a silicon — silicon-nitride — silicon waveguide 
on a silica substrate bear witness that even a three-layer stack can be designed to 
exhibit at certain frequencies the perfect phase matching of TE and TM guided modes. 
Since the TE and TM mode spectra are naturally very sensitive to refractive index 
distribution, i.e. to the index contrast of the layers and their thicknesses, number and 
succession, wide possibilities arise to tailor the mode dispersion in the vicinity of those 
phase-matching points. 

The principle of perfect phase matching with a multilayer dielectric array has 
important potential implementations in polarization-insensitive planar propagation, 
fibre-to-planar waveguide coupling, mode conversion, filtering etc. The advantage of 
the presented device scheme is its simplicity and compatibility with common integrated 
optics techniques, e.g. with planar optical waveguide technology. Further important 
issues to be addressed for the development of applicational aspects of the scheme are 
optimization with a more advanced refractive index profile of the waveguiding stack and 
integration into complex, state-of-the-art optical and optoelectronic elements. 
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Appendix 

To derive tlie mode dispersion relation for the finite, A^-period, three-layer dielectric 
superlattices, one can simply follow Yeh's ^1] reasoning for bilayer multiperiod 
waveguides, yet we prefer to utilize another tool here — the general, dyadic-based 
formalism for stratified media fH] which proves to be a competitive modification of 
Berreman 4x4 matrix approach [20] and can be directly applied to more complicated 
problems of light scattering and confined propagation in, e.g., magnetized, metallo- 
dielectric and anisotropic arrays. 

Consider an A^-period three-layer stack sandwiched between isotropic, homoge- 
neous, semi-infinite media — a substrate of permittivity esub and a cover of permittivity 
ecov The covariant dispersion relation for plane-stratified media [12] takes the form 

det (-75.6, 1) X {ViV2V^f X I ^ 

V -ICov , 

where T = z® z + y®y is the unit 2-D dyadic, z and y are the unit vectors in the plane 
of the guide, respectively parallel and perpendicular to the direction of propagation, 
a ® b denotes the outer product ajjp, Vj is the evolution operator (4x4 matrix) 
that relates z and y components of electric and magnetic fields arranged into truncated 
X X E = {—Ez, Ey)^ and Ht = {Hy, H^)^ vectors at the interfaces x = Xj and x = 
of the jth layer. 



0, (A.i; 



Hr{Xj) \p( Hr{Xj-l 

X X E{xj) J ^ \ X X E{xj^i 



-p, . . (A.2) 



and 75nfe, Icov are the impedance tensors of the ambient media that allow to generate 
given x x E: = ^ {x x E), 

'ysub,cov = 77 z(g)z^ y®y- (A. 3) 

i';Sub,Cov ^Sub,Cov 

In the case of an isotropic dielectric layer, the evolution matrix is 
= ^VjZ ® 2 + yPjy ® y 

,/% 'f^]z^zJ.. "j^.^ (A.4) 

Such a simple structure of operators Vj and 'ysub,Cov permits us to decompose 
dispersion equation (Al) into two separate equations — for TE and TM modes; these 
are obtained via multiplying ^Vi, ^V2, or ^Vi, ^V2, and applying the known 
representation of A^th power of a 2 x 2 unimodular matrix in terms of the Chebychev 
polinomials of the second kind given by Abeles |21j. Finally, this cumbersome, but 
physically transparent procedure results in equations (16)-(19). 
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